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We present a compact and simplified proof of a generalized Wick theorem to calculate the Green's 
function of bosonic and fermionic systems in an arbitrary initial state. It is shown that the decompo- 
sition of the non-interacting n-particle Green's function is equivalent to solving a boundary problem 
for the Martin-Schwinger hierarchy; for non-correlated initial states a one-line proof of the standard 
Wick theorem is given. Our result leads to new self-energy diagrams and an elegant relation with 
those of the imaginary-time formalism is derived. The theorem is easy to use and can be combined 
with any ground-state numerical technique to calculate time-dependent properties. 

PACS numbers: 05.30.-d,71.10.-w,05.70.Ln 



I 

o 



(N 

> 

00 
O 



X 



I. INTRODUCTION 

The theory of Green's functions is probably the most 
powerful and versatile formalism in physics. Due to its 
generality it has found widespread applications in any 
branch of physics dealing with many-particle systems 
such as in nuclear physics, condensed matter physics and 
atomic and molecular physics. The key of its success 
lies in the possibility of expanding the dressed (interact- 
ing) n-particle Green's function G„ in terms of the bare 
(non-interacting) Green's functions gm, m > n, which 
are then reduced to an (anti)symmetrized product of gi 
by means of the Wick theorem. -'^ In his memorable talk in 
1948 Feynman showed how to represent the cumbersome 
Wick expansion in terms of physical insightful diagrams, 
and since then the Feynman diagrams became an invalu- 
able tool in many areas of physics. 

The standard Green's function formalism (GFF) is, 
however, applicable to non-degenerate systems initially 
in their ground state and cannot be straightforwardly 
applied to systems with more general initial states or to 
systems with degenerate ground states. This leaves out 
modern fields of research like, e.g., non-equilibrium phase 
transitions,^ relaxation dynamics of ultracold gases, re- 
sponse of nanoscale systems in non-equilibrium steady- 
states,'* optimal control theory,^ etc. The correct de- 
scription of initial correlations is obviously crucial for the 
short-time dynamics of general quantum systems such as 
in transient dynamics in quantum transport^'^ or in the 
study of atoms and molecules in external laser fields. ^'^ 
In the case of finite systems it is clear that initial cor- 
relations also affect the long-time dynamics due to the 
presence of discrete quantum numbers. However, also 
in extended systems initial correlations can have infinite 
memory when the underlying Hamiltonian is integrable. 
This has, for instance, recently been demonstrated for 
interaction quenches^^"'^'''"'^^ and quantum transport^^ in 
Luttinger liquids. Another case in which the standard 
GFF is problematic, even for equilibrium properties, is 
that of systems with degenerate ground states such as in 



open shell atoms or molecules. These systems are ubiq- 
uitous in quantum chemistry but the standard Wick the- 
orem can only deal with the highest spin- or angular mo- 
mentum component of a multiplet.^'' There is therefore 
a clear need to go beyond the standard GFF. 

In 1975 Hall^* used the textbook expansion of time- 
ordered products of operators into normal-ordered prod- 
ucts and contractions to extend the Wick theorem to 
arbitrary initial states. Although the general structure is 
outlined there, the calculations of the various prefactors 
is very laborious as it requires the explicit determination 
of the sign of all possible permutations of the 2m indices 
of the gm^s. Further progress was made in the classic 
review of Danielewicz*^ in which it was shown that one 
can deal with arbitrary initial states by introduction of an 
extended Keldysh contour with an additional imaginary 
track on which the initial density matrix is represented as 
the exponential of an, in general, rt-body operator. This 
was worked out further in detail by Wagner^'^ who de- 
rived the Feynman rules for the Green's function involv- 
ing diagrams with arbitrary n-body correlators. Later 
the approach was used by Mozorov and Ropke^^ to derive 
quantum kinetic equations with arbitrary initial correla- 
tions. Furthermore Bonitz and co-workers have presented 
alternative derivations of the equations of motion of the 
non-equilibrium Green's functions with arbitrary initial 
correlations using functional derivative techniques and 
applied them to study the decay of initial correlations in 
an electron gas.^^ Recent work has further explored how 
an initial density matrix on the extended Keldysh con- 
tour can be related to diagrammatic expansions on the 
standard Keldysh contour. ^"^ 

The aim of the present work is to extend these de- 
velopments in three different ways. First of all, we rec- 
ognize that the generalized Wick theorem is simply the 
solution of the coupled system of differential equations 
for the 5m 's, the so called Martin-Schwinger hierarchy 
(MSH),^^ with proper boundary conditions. This com- 
pletely avoids the introduction of normal-ordered prod- 
ucts or contractions and therefore greatly simplifies the 
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mathematics. In particular the standard Wick theorem 
follows as a onc-linc proof. Secondly, our rcformiilation 
based on an initial value problem allows us to prove that 
the generalized Wick expansion has a form identical to 
that of a Laplace expansion for permanents/determinants 
(for bosons/fermions). Consequently, the calculation of 
the various prefactors is both explict and greatly simpli- 
fied. We provide a systematic way to express the Qm in 
terms of gi and of the initial /c-particle density matrices, 
k < m. The latter encode all the necessary information 
on the initial state and are completely determined by it. 
When the initial state is a permanent/determinant of sin- 
gle particle states the gm is expressed solely in terms of 
the gi (standard Wick theorem). Novel terms, instead, 
appear for initially correlated or entangled states. Such 
terms have a simple diagrammatic representation and 
are easier to evaluate if compared to standard diagrams. 
Thirdly, we discuss the relation between GFF based on 
the generalized Wick theorem and GFF as formulated 
on the extended Keldysh contour. These two approaches 
are complementary: the former requires the knowledge 
of the initial state while the latter requires the knowl- 
edge of the Hamiltonian with such initial state as the 
ground state. Clearly the convenience of using one ap- 
proach or the other depends on the information at hand. 
We present a straightforward proof of the existence of a 
Dyson equation for general initial states and derive an 
exact mathematical relation between its self-energy and 
the self-energy of GFF on the extended Keldysh contour. 



II. THE STANDARD GREEN'S FUNCTION 
FORMALISM 

Let us briefly review under which circumstances the 
Wick theorem can be used to calculate the n-particle 
Green's function G„. Consider a system in the state 
at time to and evolving according to the time-dependent 
Schrodinger equation with Hamiltonian H{t). At this 
stage the state is a general many-body state and 
does not have to be the ground state of H{to). We write 

H{t) = Ho{t) + W{t), (1) 

where 

Ho{t)= j dxij\x)h{x,t)i){x) (2) 

is quadratic in the field operators ip{x) and ^l^^x) (here 
a; is a one-body quantum number like, e.g., the position- 
spin coordinate) ; this decomposition is arbitrary as Ho (t) 
can be any one-body operator. The remaining term W{t) 
will often be a two-body operator describing interactions 
between the particles, but can in fact be a general n-body 
operator with any time-dependence. The basic quantities 



of the GFF are the n-th particle Green's functions 
G„(l...n; l'...n') 

= Ic^ir {v'if(i) . ..4>H{n)i'U^') . . . i^U'^')} m, 

(3) 

where 7" is the time-ordering operator, the subscript "iJ" 
denotes the Heisenberg picture and the short-hand no- 
tation 1 = {xi,ti), 1' = {x[,t'i) etc. has been intro- 
duced. The n-body Green's function gives a complete 
description of n-body correlations in a many-particle sys- 
tem. If we write the evolution operator as U{t,to) = 
Uo{t,to)F{t,to), where Uq is the non-interacting evolu- 
tion operator, then F is the evolution operator in the 
interaction picture and fulfills 

i^/{t,t') = m{t)F{t,t') (4) 

with boundary condition F{t, t) = 1 (the subscript "/" 

denotes the interaction picture). Consequently, we can 
rewrite G„ (omitting its arguments) as 

G„ = j^{^\F{to,oo)T{e-'f-'»'''^'^^M^)---i'i{n) 
X ^,\{n')...^,\{l')}F{-<x>,to)\^) (5) 

The Wick theorem applies to the g-m^s, i.e., the time- 
ordered product of field operators in the interaction pic- 
ture averaged over a non-interacting state |"I>o)- Thus to 
express G„ in terms of the gm's Eq. (5) needs to be fur- 
ther manipulated. Let us choose |$o) as the ground state 
of Ho = Ho{to) and construct the Hamiltonian 

H^ = Ho + (6) 

with ground state IvP).^^ According to the Gell-Mann and 
Low theorem^^'^^ the state \'^) (for a more mathemati- 
cal discussion on the precise conditions see Refs. 28,29) 
can be reached from |$o) by adiabatically turning on the 
interaction W^. Then, 

|*)=^^(io,-oo)|$o), (7) 

where F^ fulfills equation (4) with W p-'^lt-tol 
and 77 an infinitesimal positive constant. Equation (7) 
also implies that we can reach the state l^*) (up to a 
phase factor^^'^'') with a backward propagation in time 
and hence 

|*)=e^«F^(io,oo)|$o), (8) 

Exploiting these results the Green's function G„ in equa- 
tion (5) takes the form 

X i^|(n')...^l(l')}l^o>. (9) 
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FIG. 1: The original Keldysh contour. The contour time 
z — t± lies on the ± branch at a distance t from the origin. 
The end-points are ti = to- and tf = to+. 



where 



(10) 



In most textbooks the interaction W{t) = W is time- 
independent and [\[') = l^'o) is the ground state of 
H = H{to); then = |$o) since = W. In 

this case we can expand the exponent in equation (9) 
in powers of Wi, express G„ in terms of the noninteract- 
ing m-particle Green's functions gm's with m > n and in 
turn the g,nS in terms of 171 using the Wick theorem.'^" To 
the contrary if the interaction is time-dependent and/or 
4'} 7^ l^'o), the standard Wick theorem is of no use to 
calculate G„. Our aim is therefore to lift these restric- 
tions in a generalized version of Wick's theorem which 
can be applied to general initial states for both equilib- 
rium and non-equilibrium systems. 



III. GENERALIZED WICK THEOREM 

Reading the time-arguments in equation (5) from right 
to left it is natural to design the Keldysh contour^^ 7 of 
Fig. 1 and define the generalized Green's function 

X v-linO-.-v^la')}!*), (11) 

with %/ the contour-ordering operator and 1 — (xi,zi), 
1' = {x'i,z[) etc. collective indices with times z^s on the 
contour. The new G„ coincides with the time-ordered 
one of Eq. (5) when all contour-times lie on the upper 
branch. Expanding the exponent in powers of Wj, the 
Keldysh G„ can be written as the sum of integrals over 
non-interacting Keldysh Green's functions 

5„(l...m;l'...TO') 

= {V37(l) . . . i^i{m)^\{m') . . . 7^1(1')} I*) 

(12) 

where m > n. For these functions a generalized Wick 
theorem will now be derived. 

From the equations of motion of the field operators ipi 
and ip\ it follows immediately that the in Eq. (12) 



fulfill the non-interacting MSH 

ni 

[id,, ~ h{k)]g„, = ^(±)^+J"<5(fc,j')5rn-i(fc; j'), 

rn 

h*9.^ - h{J')]g.^ - 5](±)'=+^'<5(fc,j')ff™-i(^; j'), (13) 
fe=i 

where d{k,j') — d{xk,x'j)S{zk, z'j), h is the one-body 
Hamiltonian of Eq. (2) and the upper/lower sign refers 
to bosons/fermions. Here and in the following the sym- 
bol " " " over the indices specifies the missing indices of 
gm- Thus for instance 53(1; 2') = (73(28; 1'3'). It is worth 
noting that Eqs. (13) constitute a system of recursive 
relations since g^ can be calculated from the sole knowl- 
edge of gm-i- The solution for gm depends, of course, on 
the boundary conditions that we impose. The boundary 
conditions for gm follow directly from its definition (12) 
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(±i)™ lim .g,„(l...m;l'...m') 

2fc ,2:' -^ti 



with zi < 
hni^,_ ^ 



(14) 

< Zm < z'm < . . . < z'l (from now on 
will always be taken in such order) where ti 
is the initial time of Fig. 1 and 

^mi^\x[)...^\x'mmXm)...^{xim 

(15) 

is the initial m-body density matrix describing the ini- 
tial TO-body correlations. One can readily verify by 
expanding along a row or a column that the perma- 
nent /determinant 



5i„(l...m;l'...m') 



g{\-V) ... g(l;m') 



g{m;V) ... g{m;m') 
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of one-particle Green's functions g = gi is a solution 
of Eq. (13). However, in general the solution (16) 
will not satisfy the boundary conditions (14). This will 
only happen when the m-particle density matrix Tm of 
equation (15) is averaged over a non-interacting state 
I'?) = |<&o)- Note that this result constitutes a one-line 
proof of the standard Wick theorem which, in our formu- 
lation, amounts to solving a boundary problem for the 
MSH. 

For arbitrary initial states j^*) the particular solution 
(16) must be supplied with the additional solution gm of 
the homogeneous equations 



[id,, - h{k)]~gm = 



(17) 



and its adjoint to satisfy the boundary conditions. For 
illustrative purposes we first consider the examples of 52 



and gs- The most general solution of Eqs. (13) for m = 2 
reads 



92 = \g\2+g2- 



(18) 



Being a homogeneous solution, 52 = 52 — Iff 1 2 is not dis- 
continuous when a contour-time Zk passes through z'j, 
which implies that the discontinuity of 32 is compensated 
by an identical discontinuity in \g\2. The initial value of 
g2 is fixed by the boundary conditions (14) for (72 and g 
and reads 

lim 52(12; 1'2') - (T«)'C2(a;iX2;a;'i4) (19) 



with the correlation function 

C2 = Tz - |r|2 



(20) 



and r = Fi the one-particle density matrix. The limit 
in Eq. (19) is now independent of the time-ordering in 
which the limit is taken. The same applies to all func- 
tions (jm defined by Eq. (17), i.e., their limit is defined 
independent of the time-ordering taken. We can express 
g2 in terms of g and C2 as follows. The spectral function 
on the contour 

A{1;1') ^ *[5>(l;l')-.g<(l;l')] 

= (vl/|^,(l)^l(l')±^l(l')^/(l)l*) (21) 
satisfies the homogeneous equations 
[id,, - hil)]A{l; 1') = - hil')]Aih 1') = 0, (22) 

(hence it has no discontinuity when zi passes through z[ ) 
and equals S{xi, X2) for zi = z[. If we define 



d^{z) = S{z,ti) - S{z,tf), 



(23) 



where ti and t / are the end-points of the Keldysh contour 
of Fig. 1, we can write 



A{l,x[ti) = i[g{l,x[t,)-g{l,x[tf)] 



— i dz g{l,x'iz)S-{z) 



(24) 



and similarly 



32 



A{xiU,l') = -i / dzS_{z)g{xiz,l'). (25) 

Therefore, introducing the time-local correlation function 
on the contour 

C2(12;1'2') = S-izi)6-{z2) 

X C2{xiX2]x[x'2)6-{z[)6^{z'2) (26) 

it is immediate to realize that 

32(12; 1'2') = (Tif I dld2dVd2'g{l-l)g{2-2) 

X C2(12;f'2"')g(l';l')5(2';2') (27) 
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FIG. 2: Diagrammatic representation of g2 (a) and (b). 



provides an explicit solution to the initial value problem 
for (72, i-e., it satisfies Eqs. (17) and (19). Here and in 
the following J^dl = J dxi J^dzi. The diagrammatic 
representation of the generalized Wick theorem for 52 is 
displayed in Fig. 2(a); it is worth noting that for a non- 
interacting initial state C2 = and the standard Wick 
theorem is recovered. 

Inserting (72 = \g\2 + 92 into Eq. (13) we find in a 
similar way that the general solution for (73 is 

3 

93 = I.9I3 + E (±)'+^5(fc;j')52(fc;j') +53, (28) 

where the homogeneous solution ^3 is fixed by the bound- 
ary conditions for 53, 52 and g: 



lim 33(123; 1'2'3') = (T0't^3(a;ia;2a;3; x'l^^), 

Zk,Z -—>-ti 



with 



(29) 



3 

C3 = r3 - {±)''^'^ixk;x'^)C2{xk;x'j) - |F|3. (30) 



Using the (anti)-symmetry of 52 and 53 one can readily 
check that 53 in Eq. (28) is properly (anti)-symmetric in 
the primed and the unprimed variables. The second term 
in Eq. (28) has exactly the same structure of a general- 
ized Laplace expansion for permanents/determinants. As 
we shall see this is a general feature of all terms of the 



expansion of gm- Similarly to the two-particle case we 
can define the time-local density matrix function on the 
contour and write 

53(123; 1'2'3') = (t*)' / l[dkg{k;k) 

■^1 k=i 

3 

X C^{m-V2'i')\{dj'gQ'-j'). 

The diagrammatic representation of the generalized Wick 
theorem for g^ is displayed in Fig. 2(b). For instance, 
the second diagram of the last line in Fig. 2(b) gives 
the contribution {±f+'^ g{i,2')g2{l2]l'Z') to Eq. (28). 
Note that the sign of a diagram is given by (zb)"'^ where 
ric is the number of crossings of Green's function lines. 
Alternatively it is given by (±)"''= where ric is the number 
of transpositions needed to reorder the labels kk and j' j' 
in Eq. (28) to 123 and 1'2'3', as is readily checked for 
our example. 

At this point it is not difficult to guess the solution for 
gm. Let us introduce the collective ordered indices K — 
(ki . . . ki), with fci < . . . < fc; and Xk = {xk^ . . . XfeJ, 
and define 



Ci{K;J') = ( n S-{^kJ ] Q{Xk;X'j) I n <5-(4^ 

(31) 

where Ci{Xk',Xj) is defined by the recursive relations 

1-2 

Ci{Xk;X'j) = Ti{XK;X'j)-Y^Y.^±r^<^\ 

n=l PQ 

X \^\n{Xp]XQ)Cn-l{Xp]XQ) 

- \T\i{Xk-.X'j) (32) 

with 6*2 = r2 — |r|2. In equation (32) the sum runs over 
all ordered sequences Q = (gi . . . g„), P = {pi . . .pn) with 
indices between 1 and I, and the sign is determined by 



\P + Q\ = Y,iio 



(33) 



For example, let Z = 5 and n = 2 with P = (1,4) and 
Q = (2,4) and hence P = (2,3,5) and Q = (1,3,5). 
We then have Xp = [xi^Xi), X'q ~ {x2,x'^) and the 

complementary collective coordinates Xp — {x2,X3,xq), 
X'q = {x'i,x'^,x'^). The corresponding term for C5 under 
the summation sign in (32) is given by 



(±) 



1+4+2+4 



F 1 2 (a; 1 a;4 ; X2 0:4 ) C3 (a;2 2:3 ; x'l ) 



= ± 



r(a;i;a;2) r(a;i;a;4) 
r(x4;a;2) r(x4;a;4) 



C3{X2,X3,X5;X[,X'3,X'^). 



(34) 




+ 




FIG. 3: Diagrammatic representation of the terms of the 
expansion (35) for m — 5 and I — 2 with P = (1,4) and 
Q = (2,4) corresponding to Eq. (34). 



- Theorem: The solution of the MSH of Eqs. (13) with 
the boundary conditions (14) is 



m-2 

gm = \g\m + X! 5m +• 



(35) 



1=1 



where we defined 



.9(^(if; J') ^5](±)i^+«i|5K^;Q')3m-KAQ') (36) 

PQ 

and, in complete analogy with 32 and 33, 

I 
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{K;J') = [Ti]' f X{dk^g{k^-X) 

"'t q=1 



Ci{K;J')\{dJpg{j'p-j'p) 

/3=1 



(37) 



We will now prove the generalized Wick theorem: 



Equations (35)- (37) give a complete solution of the ini- 
tial value problem for the MSH for both equilibrium and 
nonequilibrium systems. The only input is the initial m- 
body correlations incorporated in F^. Let us make some 
additional remarks. The function gm is the sum of all 
possible diagrams with a single correlation block Cm-i in 
which m — I Green's function lines enter and leave, mul- 
tiplied by / separate Green's function lines. For instance, 
for rn = 5 and I = 2 the term of the sum in Eq. (35) 
with P = (1,4) and Q = (2,4) has the diagrammatic 
representation in Fig. 3 and corresponds to Eq. (34). 
The sign of the diagrams is (±)"^ where Uc equals the 
number of crossings (in Fig. 3 ric = 3 for the first dia- 
gram and TT-c = 6 for the second diagram). Although the 
general structure of gm was outlined in Refs. 18,19 our 
theorem contains a precise statement on how to construct 
it. When is a non-interacting product state the re- 
cursive relations in Eq. (32) yield Ci = 0, hence gi = 
and consequently gm = \g\m, i-e., we recover the stan- 
dard Wick theorem, 33 in agreement with the previous 
discussion. The generalized Wick theorem superseeds all 
common limitations of the standard GFF, i.e., the non- 
degeneracy of the ground state, the assumptions of the 
Gell-Mann and Low theorem and the requirement of a 
time-independent Hamiltonian for times t — > ±00. Fi- 
nally we would like to point out that Eq. (35) has the 
structure of the permanent/determinant of the sum of 



two matrices A and B. Indeed according to the Laplace 
formula 

m — l 

\A + BU = |A|„,+ ^5^(±)I^+Q' 

1=1 PQ 

X \A\i{P;Q)\B\rn-i{P;Q)+\B\r„ (38) 

where \A\i{P;Q) is the permanent/determinant of the 
I X I matrix obtained with the rows P and the columns Q 
of the matrix A. In the special case I = m we have P = 
Q= (l...m) and hence \A\jniP',Q) = \A\m- The same 
notation has been used for the matrix B. With the iden- 
tification Akj = g{k;j') and |B|„i_;(P; Q) = g„i-i{P; Q') 
for / = 1 ... m — 2 and the definition gi = {) Eqs. (35) 
and (38) become identical. We can thus symbolically 
write the generalized Wick theorem as 

gm = \g + g\m (39) 

whose precise meaning is given by Eq. (35). 

- Proof. The proof consists of two steps. We first 
prove that satisfies the MSH. For this purpose it 
is enough to prove the MSH for each gm since, as al- 
ready noticed, the MSH is obviously satisfied by |g|m- 
We will consider only the first of equations (13) as the 
same reasoning applies to the second. If we act with 

{idzu ~ ^(^)) '311 ^ particular diagram of gm the result 
is zero if k is an ingoing line to a C;-block. This follows 
immediately from the definition (37) and Eq. (22). To 
the contrary, if the diagram contains a free line g{k;j') 

the result is (±)'''+^(5(fc, j')x [a diagram for gm^i(k; j')]. 
Thus, {idz^. — /i(fc))^[all diagrams in gm containing 
g{k;j')] = (±)^-+^5(fc,/)g,^Li(fc;j')- The MSH for g^J} 
follows from this relation when summing over all j. Tak- 
ing into account that {idz^ — h{k))g„i = we conclude 
that gm satisfies the MSH. 

Next we prove that g^ has the correct boundary con- 
ditions (14). In Eq. (35) we take the limit Zk,Zj ti. 
Since l^]; -> (Ti)'|r|i and tjm^i {Ti)'^^^Cm-i we find 

■m-2 



Zk 




1=1 PQ 



(40) 

Writing the last term Cm as in equation (32) we see by 
inspection that all terms of the sum cancel out and the 
r.h.s. reduces to Tm- We therefore exactly satisfy the 
boundary condition (14). This concludes the proof ■ 

IV. GENERALIZED GREEN'S FUNCTION 
FORMALISM 

With our generalized Wick theorem we can now di- 
rectly expand the one-particle Green's function G = Gi 
in powers of the interaction Wj, see Eq. (11). In what fol- 
lows we restrict the analysis to interactions which are rep- 
resented by two-body operators like, e.g., the Coulomb 
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FIG. 4: First-order expansion of G for a two-body interaction. 
The last three diagrams vanish since both lines enter a Ci- 
block and hence the internal time-integration can be reduced 
to a point. This is a completely general feature. 



interaction. From Eq. (11) the interacting Green's func- 
tion is then given by 

- T.uO" / w{l,l')...w{k,k') 

X <?2fe+i(a, 1, 1' . . . fc, fc'; 6, 1+, 1'+ . . . fc+, fc'+) 

(41) 

in which w(j,j') is the two-body interaction and we 
integrate over all space-time coordinates j and j' for 
j — 1 . . .k. By inserting Eq. (35) into Eq. (41) we ob- 
tain a diagrammatic expansion of the interacting Green's 
function. For example, to first order in the interaction we 
only need the 173 of Fig. 2. If we then take into account 
that the disconnected pieces vanish we find that to first 
order in the interaction the Green's function is given by 
the connected diagrams displayed in Fig. 4. In general 
the expansion of G in powers of Wj leads to a diagram- 
matic series which starts end ends with a g-line. The ker- 
nel of this expansion is therefore the reducible self-energy 

(r) 

which we denote by Stot- the case of general initial 

(r) 

states the Sj^^-diagrams contain at most one correlation 
block and either begin and end with an interaction line 
(we call their sum !]('')) or begin/end with an interaction 
line and end/begin with a correlation block (we call their 
sum left/right reducible self-energy S(^"^^).i9 This means 
that the general structure of the diagrammatic expansion 
is 

G(l; 2) = g(l; 2) + f 5(1; i)E(^|(l; 2)5(2; 2) (42) 

where = + Y.'f^ + E^^"' and34 

J:^[\1;2)^J:'^1\1;x,)5^{z,), (43) 

eW(1;2) = <5_(zi)S«(xi;2). (44) 

Thus the self-energy is modified by the addition of a 
left/right reducible self-energy which is non-zero only if 
its right/left time-argument is one of the end-points of 



7 



+ branch 



oo 



- branch 



FIG. 5: Extended Keldysh contour with an imaginary track 
between to and to — ip. The contour time z = t± lies on 
the ± branch at a distance t from the origin and z = to — ir 
lies on the imaginary track. The end-points are ti — to- and 
tf = to - i/3. 



the Keldysh contour. The diagrammatic expansion of 
these new self-energies is very similar to that of the stan- 
dard GFF, as it is exemplified in Fig. 4. The only extra 



J 



ingredient is the appearance of the C„ 
scribe the initial m-body correlations. 



-blocks which de- 



V. RELATION TO THE IMAGINARY TIME 
FORMALISM AND EXISTENCE OF A DYSON 
EQUATION 

We next discuss the relation between the above re- 
sults and the imaginary-time formalism^^'^'' and in do- 
ing so will also generalize some recent work by Garny 
and Miiller.^'^ The average over the state \'^) can also 
be computed by tracing the field operators with p{l3) = 
Q-^H y'Tr[e~''^ ] where /? is a real positive number (in- 
verse temperature) and then letting /3 ^ oo. We re- 
mind the reader that is the Hamiltonian with \'^) as 
ground state, see Eq. (6), and can be a general n-body 
operator. This is discussed in more detail in Refs. 19-21. 
Then, 



G„(l...n;l'...n') = 



1 Tr[e-l^"°ryAt 



r dzwiiz) 



V'/(l)...^/(n)^l(n')...^|(l')}] 



Tr[e 



(45) 



r 



where je is the extended Keldysh contour of Fig. 5 and 
W{z) — for z on the imaginary track of the extended 
These new Green's functions coincide with 



contour 



35,36 



the old ones when all contour-times lie on the original 
Keldysh contour 7. The advantage of this formulation 
is that the corresponding new non-interacting Green's 
functions can be expanded using the standard Wick the- 
orem, 50, m = |5o|m, since both 50, m and go obey the 
Kubo-Martin-Schwinger boundary conditions, i.e., 
they are periodic/antiperiodic in the contour-time argu- 
ments. Note that go,m 7^ g-m- the former is averaged with 
po(/3) = e-^^VTr[e-'^^o] while the latter with p(/3). The 
new Green's function with both arguments on 7 obeys the 
equations of motion^' ^^''^^■^^ (integral over barred vari- 
ables is understood) 

[id,, ~ h{l)]G{l,2) = 5(1,2)+ / Ee(l;l)G(l;2) 



- z rsl(l;l)Gr(l;2), (46) 
Jo 

and its adjoint, where Eg is the standard self-energy, i.e., 
the same as in Wagner work,^*^ with internal integrals 
over all space and contour-times on 7e. In Eq. (46) the 
symbol "]" specifies that the first contour argument lies 
on 7 and the second on the imaginary track; the opposite 
is specified by "p'. 35^36 .^jjj ^^^^ derive an elegant 

(r) 

relation between Eg and the reducible self-energy E^^^, 
relation which will be used to prove the existence of a 



Dyson equation on 7 for general initial states. From the 
Langreth rules on 7e^^ we have 

G^il;2) = -i j dxG^{l;xto)G^{xto;2) 

+ [G^*Er.G^](l;2). (47) 

In this equation G^/^(l;2) = ±e{±ti +t2)[G>(l;2) - 
G<(1; 2)] is the retarded/advanced Green's function, G^ 
is the Matsubara Green's function with both arguments 
on the imaginary track, and we used the short-hand no- 
tation "*" for convolutions between to and to — i/3 and "•" 
for convolutions between to and 00. Inserting this result 
into the last term of Eq. (46) and taking into account 
that for any function f{xt±) = f{xt) it holds 



/(1)G(1;2) 



dxdtf{xt)G^{xt,2), (48) 



we find 



where 



[id,, ~ hil)]Gil;2) = 6il;2) 

Ee + El*G^*Er + Ee,L) -G 



(1;2), (49) 



Ee,L(l;2) = -i[J:UG^]{l;x2to)S-{z2) 
Similarly for the adjoint equation we have 



(50) 



td,,-h{2)]G{l;2)^d{l-2) 
G- fEe + El*G^*El + E 



e,R 



(i;2), (51) 
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with 



Ee,ij(l; 2) = 5-{zi)i[G'^ * ^[Kxih; 2). (52) 



Given the self-energy and the Green's function with ar- 
guments on the imaginary track wc can regard Eqs. (49) 
and (51) as the equations of motion for G with arguments 
on 7. To integrate these equations we cannot use the non- 
interacting Green's function go which satisfies the Kubo- 
Martin-Schwinger relations as the point to — ij3 does not 
belong to 7. However we can use the non-interacting 
Green's function g of Eq. (12) since it satisfies 

lim 5(1; 2)= lim G(l; 2) = r(a;i; X2). (53) 

Thus if we define the total self-energy as 

Stot = Se + Sl*G^*Er + Ee,L + Ee,K (54) 

we can write the following Dyson equation on 7 (integral 
over barred variables is understood) 



C?(l;2)=5(l;2) + 



/ 5(1; i)s 



:tot(i;2)C?(2,2). (55) 



Comparing now this equation with (42) we deduce an 
important relation between the reducible self-energy E^^^ 
written in terms of correlation blocks and the irreducible 
self-energy Etot written in terms of integral over the 
imaginary track 



EtotffEtot + EtotffEtotffEt 



(56) 



This relation constitute a bridge between two different 
and complementary approaches [each built to optimize 
the nature (density matrix or Hamiltonian) of the initial 
information] and give a deep insight in the physics of 
initial correlations. Indeed we can now express Etot in 

(r) 

terms of to obtain an irreducible self-energy in terms 
of correlation blocks. From Eq. (56) we have 



and hence 



^tot — ^t 



E^^J(l-,gEtot) 



1 



■1-ffEt, 



(57) 



= E. 



tot 



E^^SM + Ei^i^Ei^i^Si^i - . . . (58) 



Inserting this expansion into Eq. (55) we obtain a Dyson 
equation in which the irreducible self-energy depends 
only on g and correlation blocks. 



VI. SOME APPLICATIONS AND 
CONCLUSIONS 



The rules to expand the Green's function around any 
initial state constitute a powerful theoretical tool and the 

range of applicability can be vaster than we can imagine 
at present. Possible applications are in quantum chem- 
istry where spin-degenerate ground-states are ubiquitous 
while the standard Wick theorem can only cope with the 
highest spin component. The generalized Wick theorem 
is not limited to determinantal states thereby permitting 
to resolve the; spin-multiplet properties. This versatility 
can be used in the context of quantum entanglement as 
well. Another application is in the formalism of the rate 
equations which gives rise to negative time-dependent 
probabilities when the standard Wick theorem is applied 
to initially correlated states. Corrections with the C;- 
blocks should result in an improved theory. The gener- 
alized Wick theorem also allows us to address the effect 
of integrability-breaking perturbations using the exact 
Green's function of the integrable model. More gener- 
ally the generalized Wick theorem can be combined with 
any ground-state numerical technique to calculate time- 
dependent properties. Indeed from the ground-state den- 
sity matrix we can calculate the non-equilibrium Green's 
function to any order in the interaction strength. 

The practical implementation of the generalized Wick 
theorem is no more difficult than the standard Wick the- 
orem. We expect that the generalized Wick theorem will 
be useful in a broad range of physical problems both in 
and out of equilibrium. We wish to conclude by observ- 
ing that even though all derivations in this work pertain 
to Green's functions defined as averages over an initial 
state the generalization to Green's function defined as 
ensemble averages is straightforward and no extra com- 
plications arise. 
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